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We investigate the heavy quark mass effects on the parton distribution functions in the unpolarized 
virtual photon up to the next-to-leading order in QCD. Our formalism is based on the QCD-improved 
parton model described by the DGLAP evolution equation as well as on the operator product 
O l' expansion supplemented by the mass-independent renormalization group method. We evaluate 

D ( the various components of the parton distributions inside the virtual photon with the massive 

4h ■ quark effects, which are included through the initial condition for the heavy quark distributions, or 

equivalently from the matrix element of the heavy quark operators. We discuss some features of our 
results for the heavy quark effects and their factorization-scheme dependence. 
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I. INTRODUCTION 



The CERN Large Hadron Collider (LHC) has restarted Q and discoveries of signals for the new physics beyond 
the standard model (SM) are much anticipated. In order to fully complement the discoveries at LHC, we need more 
\ precise measurements which will be provided by the International Linear Collider (ILC), a proposed e + e~ collider 
machine 0. In analyzing the signals for the new physics, it is still important for us to have a detailed knowledge of 
C$ ' the SM predictions at high energies based on QCD. 

In e + e~ collision experiments at high energies, the cross section of the two-photon processes e + e~ — > e + e~ +hadrons 
dominates over other processes such as the annihilation processes e + e~ — >• 7* — > hadrons. The two-photon processes 
provide a suitable testing ground for studying the QCD predictions at high energies. We consider here the two-photon 
processes in the double-tag events, where both the outgoing e + and e~ are detected. In particular, we investigate 
the case in which one of the virtual photon is far off-shell (large Q 2 = —q 2 ), while the other is close to the mass-shell 
(small P 2 = —p 2 Y Then the process can be viewed as a deep-inelastic scattering where the target is a photon rather 
than a nucleon [3[ (see Fig. [T]). In this deep-inelastic scattering off photon targets, we can investigate the photon 
structure functions, which are the analogues of the nucleon structure functions. The study of the photon structure 
functions has long been an active field of research both theoretically and experimentally Q ■ 
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A unique and interesting feature of the photon structure functions is that, in contrast with the nucleon case, the 
target mass squared P 2 is not fixed but can take various values and that the structure functions show different 
behaviors depending on the values of P 2 . 

The unpolarized (spin-averaged) photon structure functions F^x^Q 2 ) and F2(x,Q 2 ) of the real photon (P 2 = 0) 
were studied in the parton model |5|, in perturbative QCD (pQCD) by using the operator product expansion (OPE) [6| 
supplemented by the renormalization group (RG) method @, H| and also by using the QCD improved PM Q powered 
by the parton evolution equation |U>[l3l |. The polarized photon structure function gj(x, Q 2 ) of the real photon was 
analyzed inpQCD jT34l6j . The QCD analysis has been made for F^( x, Q 2 ) u p to the next-to- next-to-leading order 
(NNLO) [Hf and for gj(x, Q 2 ) up to the next-to-leading order (NLO) 

For a virtual photon target (P 2 ^ 0), we obtain the virtual photon structure functions F^{x, Q 2 , P 2 ) and 
F2(x,Q 2 ,P 2 ). In fact, these structure functions were analyzed in pQCD for the kinematical region, 

A 2 < P 2 < Q 2 , (1) 

where A is the QCD scale parameter [l7i - l2lj ]. The advantage of studying a virtual photon target in this kinematical 
region ([I]) is that we can calculate the whole structure function, its shape and magnitude, by the perturbative 
method. This is contrasted with the case of the real photon target where in the NLO and beyond there appear 
nonperturbative pieces. With the recently calculated results of the three-loop anomalous dimensions of quark and 
gluon operators (22l . |23| and also of the three-loop photon-quark and photon-gluon splitting functions 13[ the virtual 
photon structure function F^ {F2) was investigated up to the NNLO (to the NLO) in pQCD [U,[25|. In the same 
kinematical region ([T]), the polarized virtual photon structure function gj{x, Q 2 , P 2 ) was studied in pQCD fl6l. [26l - [28| . 

In parton picture, structure functions are expressed as convolutions of coefficient functions and parton distribution 
functions (PDFs) in the target. The knowledge of these parton distributions is important since they will be used 
for predicting the cross sections of other inclusive processes. When the target is a virtual photon with P 2 being in 
the kinematical region ([1]), then a definite prediction can be made for its parton distributions inpQCD. The parton 
contents of the unpolarized and polarized virtual photon for this case were studied in Refs. fl9l l29l - [32l ] . Recently the 
QCD analysis of the parton distributions in the unpolarized virtual photon target was performed up to the NNLO [33j . 

Although the photon structure functions and photonic parton distributions were studied in pQCD up to the 
NNLO [13, HH m these analyses all the contributing quarks were assumed to be massless. The production 
channel of a heavy flavor (with mass m) opens when (p + q) 2 > Am 2 and its mass effects should be taken into account 
unless Q 2 3> m . In fact the study of heavy quark mass effects for the two- photon processes and photon structure 
functions has appeared in the literature [ill, EH, US HH l34l - [39j . Quite recently the present authors analyzed the heavy 
quark mass effects on F^ for the kinematical region ((T|) up to the NLO [4(| using a different approach from the ones 
before. The analysis was made in the framework based on the QCD improved PM and the mass-independent RG 
equations, in which the RG equation parameters, i.e., /3 and 7 functions, are the same as those for the massless quark 
case. 

In this paper we examine the heavy quark mass effects on the parton distribution functions in the unpolarized 
virtual photon up to the NLO in pQCD. We use the same framework as the one in Ref. [4(|, the QCD improved 
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PM combined with the mass- independent RG equations. We consider the system which consists of n/ — 1 light (i.e., 
massless) quarks and one heavy quark together with gluons and photons. Then, the heavy quark mass effects are 
included in the RG equation inputs; the coefficient functions and the operator matrix elements. In the case of the 
nucleon target, the heavy quark mass effects were studied by a method based on the OPE in Ref. j4l[, where the 
heavy quark was treated such that it was radiatively generated and absent in the intrinsic quark components of the 
nucleon. This picture does not hold for the case of the photon, since the heavy quark is also generated from the photon 
target together with light quarks at high energies. We should consider both the heavy and light quarks equally as the 
partonic components inside the virtual photon. 

In the next section, we derive the evolution equations for the parton distribution functions for the case where nj — 1 
light quarks and one heavy quark are present. Then solving these equations, we give the explicit expressions up to 
the NLO for the moments of the light flavor singlet (nonsinglet) quark, heavy quark and gluon distributions. The 
parton distributions are dependent on the scheme which is employed to factorize structure functions into coefficient 
functions and parton distributions. We investigate the photonic parton distributions in two factorization schemes, 
namely MS Hi and DIS 7 HH schemes. In Sec. EH we enumerate all the necessary QCD parameters to evaluate the 
photonic parton distributions up to the NLO in MS scheme. The parton distributions in DIS 7 scheme are considered 
in Sec. IIV1 The numerical analysis of the parton distributions predicted by MS and DIS 7 schemes will be given in 
Sec. El The final section is devoted to the conclusions. In Appendix [A] we consider the parton distributions in the 
virtual photon for the case when all rif quarks are light. 



II. PARTON DISTRIBUTIONS IN THE VIRTUAL PHOTON WITH A HEAVY QUARK FLAVOR 



We investigate the parton distributions in the virtual photon for the case when one heavy flavor quark appears 
together with n/ — 1 light (i.e., massless) quarks. The analysis is made in the framework of the QCD improv ed p arton 
model @ powered by the DGLAP parton evolution equations. A part of the analysis was reported in Ref. |40| . Let 



ql(x,Q 2 ,P 2 ), ql{x,Q\P 2 ), G~>(x,Q 2 ,P 2 ), r 7 (x,Q 2 ,P 2 ), 



(2) 



be light quark (with i-flavor and i = 1, • • • , rif — I), heavy quark, gluon and photon distribution functions in the 
virtual photon with mass — P 2 . Since the parton distributions in the photon are defined in the lowest order of the 
electromagnetic coupling constant, a = e 2 /47r, T 7 does not evolve with Q 2 and is set to be r 7 (x, Q 2 ,P 2 ) = 6(1 — x). 
In the light quark sector, it is more advantageous to treat, instead of using q % L , the "flavor singlet" and "nonsinglet" 



combinations 



and 



defined as follows: 



lis =J2 q L' 

i=l 
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iLns 
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rif — 1 
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<1l, 



(3) 



where is the electromagnetic charge of the i-flavor quark in the unit of proton charge. 
Now introducing a row vector 

<f = (?Z,j«h>g !7 .?L.)> 



(4) 



the parton distributions q 1 (x, Q 2 ,P 2 ) in the virtual photon satisfy the inhomogeneous DGLAP evolution equation [Tol — 



rig 7 (x,Q 2 ,P 2 
d\nQ 2 



k(x,Q 2 



l dy 



^(y,Q',P 2 )P[-,Q 



y 



y 



(5) 



where the elements of a row vector k = (k^ s , fc# , kc k^ ns ) refer to the splitting functions of 7 to the light flavor-singlet 
quark combination, to heavy quark, to gluon and to light flavor-nonsinglet combination, respectively. The 4x4 matrix 
P (z, Q 2 ) is expressed as 



P(z,Q 2 ) = 

where Pab is a splitting function of B-parton to A-parton 



(P s LL {z,Q 2 ) Phl(z,Q 2 ) Pgl(z,Q 2 ) 

Plh(z,Q 2 ) P H h(z,Q 2 ) P gh {z,Q 2 ) 

Plg(z,Q 2 ) Phg(z,Q 2 ) Pgg(z,Q 2 ) 




pNS 







(z,Q 2 )J 



(6) 
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The method to solve the above inhomogeneous DGLAP Eq. ([5]) is well known ll|. We sketch out the procedures. 
First we take the Mellin moments, 



dq~<{n, Q 2 ,P 2 ) 
d\nQ 2 



k(n, Q 2 ) + q 7 (n, Q 2 , P 2 )P{n, Q 2 ) 



(7) 



where we have defined the moments of an arbitrary function f(x) as f(n) = dxx n 1 f(x). Hereafter we omit the 

obvious n dependence for simplicity. Then expansions are made for the splitting functions k(Q 2 ) and P(Q 2 ) in powers 
of the QCD and QED coupling constants as 



HQ 2 ) - £ fc (0) + ^M 1 ^ 



P(Q 2 ) 



a s (Q 2 ) p(Q) 
2tt 



a s (Q 2 ) 



2vr 



and a new variable t is introduced as the evolution variable instead of Q 2 [44j , 

'"A a s {Q 2 ) ■ 

The solution q 1 {t)(= q 1 (n,Q 2 ,P 2 )) of ([7]) is decomposed in the following form: 

q^(t) = q^ \t) + q^(t) , 



(8) 
(9) 

(10) 
(11) 



where the first and second terms represent the solution in the LO and NLO, respectively. Then they satisfy the 
following two differential equations: 



dq~<(°\t) 


a 




dt 


a s (t) 
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dt 


2^ 





^_ k (0) 



2/?o 



2tt 



,7(0) 



P (i) _ 

2/?o 



? 7(D W p(0) 



where we have used the fact the QCD effective coupling constant a s (Q 2 ) satisfies 



da s (Q 2 ) R a s (Q 2 ) 2 R a s (Q 2 ) 3 , 
— vf~ = ~ Pi ~ ^ r 



rflnQ 2 



4tt 



(4tt) 2 



(12) 
(13) 

(14) 



with Po = (11 — §H/) and fix = (102 — ^-rtf). Note that the P 2 dependence of q 1 solely comes from the initial 
condition (or boundary condition) as we will see below. 

The initial conditions for q 1 ^ and q 1 ^ 1 are obtained as follows: for — p 2 = P 2 A 2 the photon matrix elements 
of the hadronic operators O" (i = qi, H, G, Lns) can be calculated perturbativcly. (These hadronic operators O™ are 
explained in Sec. IIII A|) . Renormalizing at a 2 = P 2 , we obtain at one-loop level 



(l(p) I 0?(n) I 7(P))U»=J« = , * = QL,H, G, Lns . 



(15) 



The An terms represent the operator mixing between the hadronic operators and photon operators in the NLO and 
the operator mixing implies that there exist parton distributions in the photon. Thus we have, at a 2 = P 2 (or at 
i = 0), 



,7(0) 



(o) = o, ^) ( o) = ^2i 1) , 



with 



(16) 



(17) 



which state that the initial quark distributions emerge not in the LO (the order a/ a s ) but in the NLO (the order a), and 
initial gluon distribution starts to emerge in the NNLO (the order aa s ). Despite of the initial condition q 7 ^°'(0) = 0, 
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the LO quark distributions, both light and heavy, are generated from the photon distribution T"'(x, Q 2 , P 2 ) = 6(1 — x) 
(see Eq. ([2])) through the pointlike coupling of the photon to quarks. The heavy quark parton appears in the LO as 
a massless quark, while, as we see in Sec. IIIIB1 the heavy quark mass effects arise from the initial condition g 7 ^ 1 ) (0) 
more closely, An (1) in Eq. (32)) in the NLO. 
With these initial conditions p^|) . the solutions q 1 ^ '^) and q 7 ^ 1 ^) are given by 



,7(0) 



,7(1) 



(t) 



(t) 



4-7T 

-2 a 
+b J 1- 



a,(0) 



dre^ 

a s {t) 
a.(0) 



(0)_/3o s 



P (l 



3o 



Ap(o) 
2/?o 



,7(1) 



(0) 



a s (t) 



a„(0) 



_p(°> T 



(18) 



(19) 



where 



2tt(3 



. k (o). 



1 



b = ^ 



a 

2tt 



_ 2P<°) ' 

k 2/3 k 



2 a 



p(i) _ Ap(«) 
2/?o 



p(o) ' 



(20) 
(21) 



The moments of the splitting functions are related to the anomalous dimensions of operators as follows: 



p(0) 



1 



9(0) 

In i 



p(i) 



l 



9(1) 

In i 



i n i 



8 ™ ' 



(22) 
(23) 



where 7^°^ and 7n are the 4x4 one-loop and two-loop anomalous dimension matrices in the hadronic sector, 
respectively, and K„ an d are the four-component row vectors which represent the mixing in one-loop and 

two-loop level, respectively, between the photon operator and the four hadronic operators. The details are explained 
in the next section. 

The evaluation of q^°\t) and <j 7(1) (i) in Eqs. (18) and (19) can be easily done by introducing the projection 
operators P" such as 



p(0) 



pn pn _ 

i 3 



: 7° 

1 in 



o 

PI 



1 
4 

i ^ 3, 
i = i 



+ , — ,Lns 



Vv 



-, Lns, 



i=?^ ) -(- ) — } Lns 



(24) 
(25) 



where A™ are the four eigenvalues of the matrix 7°. Then, rewriting a s (0) and a s (f) as a s (P ) and a s (<3 ), respectively, 
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we obtain 



,7(0) 



(*)/ 



8tt/3 



4?r 



a s (Q 2 ) 



< 0) E^ r 



i - 



a s (Q 2 



2\-\ !+<*? 



,(P 2 ) 



(26) 



,7(1) 



(*)/ 



8tt/3 



= J V P" — + ^llfW VP." (l - — 



pra C;(l) pn 



n 4f 2 /3 + A; 1 - X? d. 



_ _L _ 2^ a1 1) y pj> 



2/3oAi 1} , 



7n 



X 1- 



a s (Q 2 ) 



2/3 + A™ 



Ay i 



x 1 



a,(Q 2 ) 



a s (P 2 



where d™ = and i, j = ij), +, — , Lns. 

Finally, since q 7 (t) = q 7 ' )(i) + q 1 ^ 1 '^), and from ((4]), the moments for the parton distributions of the "flavor- 
singlet" light quark, heavy quark, gluon and "fiavor-nonsinglct" light quark arc given, respectively, by 

q1 s (n, Q 2 , P 2 ) = (1, 1) component of the row vector g 7 (i) , 
q]j(n,Q } ,F ) = (1,2) component of the row vector g 7 (i) , 
G 7 (n, Q 2 ,P 2 ) = (1,3) component of the row vector g 7 (i) , 
q2 ns (n,Q 2 , P 2 ) = (1,4) component of the row vector g 7 (i) . 



(27) 



(28) 
(29) 
(30) 
(31) 



III. PARAMETERS IN q 7(0) (£) AND q 7(1) (i) 



We give here the information on the parameters which appear in g 7 (°'(i) and g 7 *- 1 ^) in ([26]) - ([27)) . They are 
calculated in MS scheme [42|. We introduce the following quark-charge factors in the massless quark sector: 



, n f -l n/-. 

— 1 ^-^ nt — 1 ^-^ 

i=1 J 1 = 1 



n f — 1 



(32) 



A. Anomalous dimensions 



Corresponding to the splitting functions given in Eq. ([6]), the anomalous dimensions in the hadronic sector are 
expressed in the form of 4 x 4 matrix as 



7n(fl) = 



l n qLH {9) 1hh{9) 7g ff (ff) o 
7™ iG (ff) 1hg(9) 1gg(9) 



V ll ns {g)J 

The four-component row vector 

K n (g,a)=(K^(g,a), K${g,a), K%(g,a), Kl ns {g,a)) , 



(33) 



(34) 
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represents the mixing between photon and four hadronic operators. Here an importance of inclusion of the heavy 
quark operator should be stressed. We treat the heavy quark in the same way as the light quarks and assume that 
both heavy and light quarks are radiatively generated from the photon target. In contrast, in the case of the nucleon 
target, heavy quarks are treated as radiatively generated from the gluon and light quarks. 

Since the elements and 7]y gi start at the order of a 2 and J^ l q has a factor {rif — 1), the one-loop anomalous 

dimension matrix 7^ is expressed as 



;(o) 







(0),n 

(0),n (0),n 

nj-l (0),n J_A0),n (0),n 
n f ItpG n f >4>G <GG 



(0),- 

















(35) 



7W7T , 



where 7^'™, 7^3 ™ ', 7g1'™ anc ^ 7gg are well-known one- loop anomalous dimensions for the hadronic sector which 
appear when all n/ flavor quarks are massless, and are given, for example, in Eqs. (4.1), (4.2), (4.3) and (4.4) of 



Ref. 



J . The four eigenvalues of jn are 

yi _ (0),n 

A 1/l — ll/)^ ' 



A Lns 



(0),' 



(36) 
(37) 
(38) 



The one-loop anomalous dimension matrix 7^ can be expressed in terms of its eigenvalues A" (i = ip, +, — , Lns) and 
corresponding projection operators as 
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\nl 1 cJ°)'' 
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n/ 7i/jG 
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7 G ^ 

(0),n 

7gg 



An 



pn 

Lns 



/0 
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(40) 



(41) 



(42) 



The elements of the one-loop anomalous dimension row vector K) t 



(0) _ f K (0),n K (0),n K (0),n ^(0),n 



H 



G 



are given by 



K (o). 



K 



(0); 

H 



K (o). 

A G 



K 



(0),n 
Lns 



24(n f -l)(e 2 ) L k° n , 

24e 2 H k° n , 

0, 

24(n / -l)((e 4 > i -((e 2 ) i ) 3 ) *£, 



(43) 
(44) 
(45) 
(46) 



with 



n 2 + n + 2 
n(ra + l)(n + 2)' 



(47) 
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The two-loop anomalous dimensions for hadronic sector with a heavy quark are inferred from those for the case 
when all n/-flavor quarks are massless (45[. Minor changes of group factors arise from quark loops: 



(l),n 
^ Lns 



(l),n 
Tq£G 

(l),n 

(l),n 
7i?ff 

(l),n 
^HG 

(l),n 
(l),n 

~1gh 



(l),n 

7jvs 



7ws 



8^(^) fl^ + gCx^)^ 



( - ) Dp S ^, , 



S"+o(i)i> 



7jv<? 



(l),n 
7gv> 

(1),™ 
Tg^ 



PS,tl>tl> > 

1 

2 



(48) 
(49) 

(50) 

(51) 

(52) 

(53) 

(54) 

(55) 
(56) 



where Ca = 3 and C F — § in QCD. The anomalous dimensions 7^" and 7qA'" for the case of n/ massless quarks are 



given, for example, in Eq. (3.5) of Ref. [22j and Eq. (3.8) of Ref. [23||, respectively (see also Ref. j24|). The expression 
of the "pure singlet" contribution D ps ^ and those of D^ G and E^ G (which appear in the two-loop anomalous 

dimension in the case of rif massless quarks) are given in Eqs. (3.6) and (3.7) of Ref. (23j, respectively. The 

anomalous dimension 7^'" for the case with a heavy quark is the same with the case of rif massless quarks and is 



given in (3.9) of Ref. \2c 

The elements of the two-loop anomalous dimension row vector = ( ifg* , Kjp' n , K G ^' n , Kj^'J 1 ) are given by 



K$> n = -3(n f -l)(e 2 ) L C F 8D 



iiGi 



K 



K, 



(l),n 
H 

(l),n 



= -3 e 2 H C F 8ffi 



K 



where D GG is obtained from j GG n by replacing Ca — > and C F nf — > \. 



= -3 ((«/ - l)(e 2 )L + e 2 H )c F 8(D GG - 1), 
= -3(n, - 1) «e 4 ) L - «e 2 ) L ) 2 ) C F 8^ G , 



(57) 
(58) 
(59) 
(60) 



B. The one-loop photon matrix elements 



The elements of the row vector A n in Eq. (fTT)) are given by 

A*?™ = 3(n / -l)( e 2 ) LJ ffW(n), 



where 



ffW(n) 



n 2 + ?i + 2 



TLns(l) 

-5i(n) 



3(n / -l)((e 4 ) L -(( e 2 ) L ) 2 )i/( 1 'H, 



3e 



n(n + l)(n + 2) 

n 2 + n + 2 



(n + 2) 2 (n + 1) 2 



In- 



1 1 



n 2 + n + 2 



;Si(n) 



(61) 
(62) 
(63) 

(64) 
(65) 



n(n+l){n + 2) P 2 n n 2 (n + l) 2 (n + 2) 2 n(n + l)(n + 2) 

with 61(71) = y,"_ 1 j, and A^' 1 ^ is obtained by evaluating the diagrams in Fig. din the limit A 2 < P 2 < to 2 . The 
heavy quark mass effects reside in the term AA 7 ^j. 
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FIG. 2: The diagrams for An . The double lines express the heavy quark. 



IV. THE NLO PDFS IN DIS 7 SCHEME 

The structure functions of the photon (nucleon) are expressed as convolutions of coefficient functions and parton 
distributions of the target photon (nucleon). But it is well known that these coefficient functions and parton distri- 
butions are by themselves factorization-scheme dependent. The relevant quantities given in Sec. Mil were calculated 
in MS scheme. When we insert them into the formulae given by (|26| and (|27|) . we obtain the parton distributions 
predicted by MS scheme. Meanwhile, an interestin g an d also useful factorization scheme called DIS 7 was introduced 
in the analysis of the photon structure function K 7 [43|. In this scheme, the photonic coefficient function , i.e., the 
direct photon contribution to F%, is absorbed into the photonic quark distributions. 

The Mellin moments of the virtual photon structure function -F^fa, Q 2 ,P 2 ) is expressed as 

/ dxx n ~ l - F?(x,Q 2 ,P 2 ) = F^{n,Q 2 ,P 2 ) 
Jo x 

= q*{n, Q 2 , P 2 ) ■ C 2 (n, Q 2 ) + C](n, Q 2 ) , (66) 

where q 7 (n, Q 2 ,P 2 ) is the four-component row vector given in @. The column vector C 2 (n, Q 2 ) is made up of four 
hadronic coefficient functions, 

C 2 (n,Q 2 ) ee (C^(n,Q 2 ), Qf(n,Q 2 ), C 2 G (n,Q 2 ), C 2 L ">, Q 2 )) T , (67) 

where C^ s , , C§ and C 2 Lns are coefficient functions corresponding to the light "flavor-singlet" quark, heavy quark, 
gluon and light "flavor-nonsinglet" quark, respectively. The last term C^n, Q 2 ) in (|66|) is the photonic coefficient 
function. The moments of the parton distributions in DIS 7 scheme are obtained as follows [H|,|43[. In this scheme, the 
hadronic coefficient functions are the same as their counterparts in MS scheme, but the photonic coefficient function 
is absorbed into the quark distributions and thus set to zero, 

C 2 (n,Q 2 )\ Bls _ i =C 2 (n,Q 2 )\ m , C 2 >, Q 2 )|dis 7 - 0. (68) 

Then Eq. §^ gives 

F](n 7 Q 2 ,P 2 ) = g>,Q 2 ,P 2 )| D iS 7 • C 2 (n, Q 2 )| DI s 7 

= q>,Q 2 ,P 2 )| DI s 7 -C.KQ 2 )!^ . (69) 
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On the other hand, F^in, Q , P ) is expressed in MS scheme as 



F 2 >, Q 2 , P 2 ) = g*r(», Q 2 , P 2 )\ m ■ C 2 (n, Q 2 )^ + C 2 >, Q 2 ) 



MS 



(70) 



We expand <j» 7 (n, Q 2 , -P 2 )|dis t in terms of the LO and NLO distributions and Ci{n, Q 2 )\m$ m powers of a s (Q 2 ) up 
to the NLO as follows: 



q~<(n,Q 2 ,P 2 



DIS. 



^ (0) +^ (1) |dis, 



^2(n,Q Jlj^g - C 9 l , H ; C 



'2,ra 



47T 



2,nlMS 



(71) 
(72) 



where the LO qZ^ and C^' n are both factorization-scheme independent. Denoting the difference of qZ {1> \ms-, from 



,(0) 



7(1) I 



MS scheme prediction as 5qZ^ |dis t , we write 



,7(D | =,7(1)^ + ^7(1)1 



Putting dTTJ)-((73j) into the r.h.s. of ([69]) and comparing the result with (f70|) . we find 

C](n,Q 2 )\ m = SqZ^\ Dls , i -C^ + --- . 



Since 



the r.h.s. of Eq. (|74]) is rewritten as 



,(0) 



((e 2 ) L ,e 2 J ,0,l) T , 



,(0) 



,7(1) 



2 A„7(l)| 



.7(1) 



The moment of the photonic coefficient function C 7 (ra, <5 2 ) Ims * s written up to the one-loop level as 

y](n, Q 2 )\y^ 



C 2 >, Q 2 )^ = ^3{( n/ - l)( e 4 ) L i?^" + 4^"} + • • 



Now dividing the light quark-charge factor (e 4 ) l into two parts, the light flavor-singlet and nonsinglet parts, as 

(e 4 )L = (e 2 )L(e 2 )L + ((e 4 )L-((e 2 ) L ) 2 ) , 
and from Eqs. (76]) and jTTJl we find at the NLO 

^ZKlmS, = ^3(n f -l)(e 2 ) L B^ , 



A«7(l)| 

|dis 7 

.7(1) 



4tt 



3 e% B"' n , 



«i,„|Dis 7 - -3(n f -l)((e^ L ^((e 2 ) L ) 2 )B 



The coefficient B^' n is related to the one-loop gluon coefficient b'q by B^' n = ^jBg @i anc ^ gi VGn by 



B 



L.n 



n 2 + n + 2 . a i w 4 
(1 + S 1 (n)) + 



n(n + l)(n + 2) 



(n + 1) (n + 2) n 2 



while B^'™ is calculated in the heavy quark mass limit (A 2 <C P 2 -C to 2 ) and we find 



B 1 -B 1 



Finally in DIS 7 scheme we set in all orders 



G 7 (n,Q 2 ,P 2 )| DIS = G 7 (n,Q 2 ,P 2 ) 



MS 



(73) 
(74) 
(75) 

(76) 

(77) 

(78) 

(79) 
(80) 
(81) 

(82) 

(83) 
(84) 
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V. NUMERICAL ANALYSIS FOR PDFS WITH HEAVY QUARK EFFECTS 



The parton distributions are recovered from their moments by the inverse Mellin transformation. Using the formulae 
given in Eqs. (|2G1) - (|3T1) and parameters enumerated in Sec. IIII1 wc obtain the parton distribution functions in the 
virtual photon in MS scheme up to the NLO . We have considered the following two cases: 

(i) Q 2 = 5GeV 2 and P 2 = 0.35GeV 2 , 

(ii) Q 2 = 30GeV 2 and P 2 = 0.35GcV 2 . 

In both cases we take n/ = 4, and choose c as a heavy quark and assume that the other u, d and s-quarks are massless. 
We take m c = 1.3GeV as an input for the charm quark mass and put A = 0.2GeV for the QCD scale parameter. 

The values of Q 2 and P 2 for the case (i) correspond to those of the PLUTO experiment (4^|. We plot the parton 
distribution functions in MS scheme in Fig. [3J for the case (i) and in Fig. Ufor the case (ii): (a) the light "flavor- 
singlet" quark distribution xq1 s (x, Q 2 , P 2 )\- W jg, (b) the heavy (charm) quark distribution xq] 1 (x, Q 2 , P 2 )\-&Js, (c) the 
gluon distribution xG 7 (x,Q 2 ,P 2 )\^jg and (d) the light "flavor-nonsinglet" quark distribution xq] jns (x, Q 2 , P 2 )\y^. 
In order to see the heavy quark effects, we plot, in addition, the parton distributions xq1 s \ light Ms> x 9h I light Ms an< ^ 
^C 7 liight ms which are obtained when c-quark is also set to be massless. Actually, we get these distributions by 

setting AAfj — > in Eq. (|62|) and by inserting the "new" row vector A^ n ^ into the expression of Eq. (|27jl . See also 
Appendix [A] Since the light "flavor-nonsinglet" quark decouples to the other partons, we have x<72„Jiight = xq 1 Lns . 

We observe from Fig. [3J (a)-(c) and Fig. 21 (a)-(c) that the c-quark mass has rather large effects on the c-quark 
distribution xqjj |jjg and the gluon distribution xG 1 \ygr, while it has negligible effects on the light "flavor-singlet" quark 
distribution £<zZJms- The difference between x<z2Jms anc ^ x< ?Lsliight MS ( an d also between xqjj\y^ and xg^| light j^g 
and between xG 1 \y^ and £G 7 | light jjg) is due to the appearance of AA]j in Eq. (|62|) . which is negative and larger in 

magnitude for smaller n. The negative means that once c-quark has mass, the c-partons are less produced from 

the target photon. Actually it has the same effect as reducing the evolution of c-quark distribution to the range of 
to 2 to Q 2 instead of P 2 to Q 2 . Indeed we find from Eq. ((27)) . 








/ \ \ f (0),n , n (0),n * n . 

<L* ~ «l„Jw) I (j^) - (l - £) - irf> + \^(r)* + ^-^M^ } , ,85) 



L\ M<Z x 1 7 ^ I A - M< x 1 7 ^ ~ A + 



(86) 



where r = ^pjpry- Unless n is a small integer, we see 7^'™ ~ A™ and d", « d™. Therefore the sum in the curly 
brackets of Eq. diminishes, which means that the effects of heavy quark on xqJ jS are extremely small. Sec Fig. [3J 
(a) and Fig. 21 (a). On the other hand, the sum in the curly brackets of Eq. (|86p is expressed approximately as (r) d * 
for n not being a small integer. The ratio of (qjj n — qjj J light) to the leading order qjj°2 is proportional to the product 

of a s (Q 2 ) and (r) d *. The values of r and a s (Q 2 ) are 0.463 and 0.237, respectively, for the case (i) and 0.351 and 
0.180 for the case (ii). The ratio is not small and Fig. [3] (b) and Fig. [4] (b) show the large reduction of xq^-^m from 

^nliight.MS' es P eciall y> for the case (i)- 

The gluons do not couple to the photon directly and they are produced from the target photon through quarks. 

Therefore, the leading contribution to the gluon distribution xG 1 1 jjjg is essentially of order a and it is very small in 

absolute value except in the small x region. The c-quark mass effects appear in xG 1 in the NLO (the order a) and 

are enhanced by the charge factor en = 2/3 (see Figs. |3j (c) and|4] (c)). The departure of z<z2sIms lrom ^Zsliight Ms 

at small x is related to the behaviors of the gluon distributions xG 1 \y^ and ^G 7 || ight ^jg. As x — > 0, the both gluon 

distributions grow while their difference becomes larger. 

Figs. E](d) and 2](d) show the light "flavor-nonsinglet" quark distribution xq1 ns \^g. Comparing with the graphs of 
x 1~Ls ImS' it i s ver y sma ll m absolute value. This is due to the fact that xq1 ns |^jg has the charge factor ((e 4 ) i — ((e 2 ) l) 2 ) 
which is a very small number (2/81 for rif = 4). An examination of Fig. [3J (b), (c) and Fig. 21(b), (c) show that with 
larger Q 2 , the c-quark mass effects become smaller. When Q 2 gets still larger, we may need to consider the 6-quark 
mass effects with taking rif = 5, but 6-quark has milder effects than c-quark because of its charge factor. 

We see from Fig. [3J(a), (b), (d) and Fig. 2] (a), (b), (d) that the quark distributions £<zZJms"' xqjjly^ and xq1 ns \j^g 
diverge as i — > 1. This is due to the NLO contributions to the quark parton distributions in MS scheme. The 
behaviors of parton distributions near x — 1 arc governed by the large-n limit of those moments. In the leading order, 
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parton distributions are factorization-scheme independent. For large n, the moments of the LO quark distributions, 
q]}g \ q]}°^ and <?2ni' behave as l/(n Inn). Thus, in x space, these LO quark distributions vanish for x — » 1 as 
[— l/ln(l — a;)]. On the other hand, the moments of the NLO quark distributions in MS scheme, ^^ImS' Qh^Ims 
and ?2„ s ImS' behave in large- n limit as (In n)/n. Therefore, in x space, the (LO+NLO) quark distributions in MS 
scheme positively diverge as [— ln(l — a;)] for x — > 1. The moments of the LO and NLO gluon distributions, G7/°' ) and 
Gn , behave for large n as l/(n ln?i) 2 and l/n 2 , respectively, and thus, in x space, the (LO + NLO) curve of the 
gluon distribution (both G 7 |j^g and G 7 | light j^g) vanishes as (— lnx) for x — > 1. 

In DIS 7 scheme the photonic coefficient Gj is absorbed into the quark distributions. In consequence, the (LO+NLO) 
quark distributions show different behaviors at large x from those in MS scheme. Since B^> n {= B^' n ) in Eq. ([77jl 
behaves as (— 41nn)/n for large n, the (LO+NLO) curves in x space for qlJms^, <?#|dis.v and <?2„ s |dis 7 negatively 
diverge as ln(l — x) for x — Y 1. In fact, using Eqs. (|79[) - H84H and inverting the moments, we obtain parton distributions 
in DIS 7 scheme up to the NLO, which are plotted in Fig. [3] and Fig. [S] Again we have considered the two cases: 
(i) Q 2 = 5GeV 2 , P 2 = 0.35GeV 2 and (ii) Q 2 = 30GeV 2 , P 2 = 0.35GeV 2 . The other parameters are the same as 
before and c-quark is taken to be heavy. We see from Fig. [5] (a), (b), (c) and Fig. [S] (a), (b), (c) that the quark 
distributions xg2 s |Dis 7 j xQhIdis., an d £9Z ns |Dis T become negative at large x. We observe again that the mass of c- 
quark has negligible effects on the light "flavor-singlet" quark distribution xqljuis-, but large effects on the c-quark 
distribution an^lois,,- When Q 2 gets larger, the heavy quark mass effects become smaller. It is noted that if we take 
into account the charge factors, the following three "renormalized" distributions, £<z2Jii g ht,Dis T = 2 ; <?2 s |iight,Dis 7 /(e 2 )i, 
^9^|ii g ht,Dis T = £<7# liight.Dis^/e^ and xql ns \uis-, = xqlnslvis-, / ((e 4 )i - ((e 2 )i) 2 ) overlap for almost the whole x 
region except near x = (see Fig. [5] (a), (b) and (c) and Fig. H](a), (b) and (c) ). 

Finally the gluon distribution :tG 7 |dis is the same as xG^lj^g. 



VI. CONCLUSIONS 



We have studied the heavy quark mass effects on the parton (light singlet, heavy quark, gluon, light nonsinglct) 
distribution functions in the virtual photon up to the NLO in perturbativc QCD. Our calculation is based on DGLAP 
equation as well as on the OPE formalism within the framework of the mass-independent renormalization group. The 
heavy quark effect is included through the operator matrix element for heavy quark operator and is evaluated by 
the heavy quark mass limit (A 2 <C P 2 <C m 2 ). In the language of the parton picture, the heavy quark effects arc 
arising from the initial condition for the heavy quark distribution. In fact, the leading- logarithimic term of our initial 
condition (|65[) can be reproduced by solving the boundary condition, qJ^(x,Q 2 — m 2 ) = [3lL l35|. in the leading 
order approximation. 

The heavy quark mass effects tend to reduce the values of parton distribution functions for the light-singlet, the 
heavy-quark and the gluon distributions except for the light nonsinglct distribution. Especially the suppression for 
the heavy parton distribution for the up- type quark is enhanced by the dependence of the charge factor = 2/3. 
These behaviors are consistent with our previous work on the virtual photon structure functions with the heavy quark 
mass effects These results could be explained by the suppression of the evolution range due to the mass of the 
heavy quark. We have also studied the factorization-scheme dependence of our parton distributions with heavy quark 
mass effects, especially for two factorization schemes, MS and DIS 7 . 

In our formalism where we treat the contribution from the twist-2 operators to the parton distributions, we have 
not taken into account the kinematical threshold effects which manifest as the presence of the maximal values of the 
Bjorken variable. We need some improvement in which the threshold effects are included. We should also investigate 
the general kinematical region where P 2 and m 2 are of the same order. We also note that the general-mass variable- 
flavor-number scheme (GM-VFNS), which has now become a popular framework for the global analyses of parton 
distributions, should be implemented in the present analysis which is under investigation. 

Such improvements would help us to predict the photonic parton distribution functions which could be measured 
at the future linear collider ILC. Especially, the application of our results to the up-type quark parton distribution 
functions like the charm quark, the top quark in the unpolarized virtual photon will be important phenomcnologically. 
The application of our formalism to the polarized photonic parton distribution functions can be carried out and would 
turn out to be relevant for the measurement of the polarized photonic PDFs at ILC. 
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FIG. 3: Parton distributions in the photon in MS scheme for n } = 4, Q 2 = 5GeV 2 , P 2 = 0.35GeV 2 with m c = 1.3GeV 
and A = 0.2GeV: (a) xq[ s {x, Q 2 , P 2 )|ms and x <lL ( x , Q 2 , ^ 2 )liight,Ms; ( b ) x 1b{ x , Q^ p2 )\ms and x Qh( x > Q 2 > -P 2 )ln g ht,MS' ( c ) 
xCnfaCf,!*)^ and xC (x,Q 2 , P 2 )\ lightMB ; (d) xql ns (x, Q 2 , P 2 )\^. 



Appendix A: PDFs for the case of massless quarks 



In this paper wc have considered the parton distributions in the virtual photon for the case when the nf-th flavor 
quark is heavy and the rest of rif — 1 flavor quarks are light (i.e., massless) and we have derived the formulae for the 
moments of the parton distributions, q1 s (n,Q 2 ,P 2 ), q]j(n, Q 2 , P 2 ), G^in^Q 2 , P 2 ) and q1 ns {n,Q 2 , P 2 ) up to NLO, 
which are given in Eqs. (p^ - (j3"T|) . If the rc./-th flavor quark is also light, in other words, all the n/ flavor quarks are 
light, we obtain, instead, the parton distributions of q1 s (n 7 Q 2 7 P 2 )|n g ht, Qnin, Q 2 , -P 2 )|n g ht and G 7 (n, Q 2 , P 2 )|ii g ht 
and q~L ns (n, Q 2 , -P 2 ) liight- Here it is noted that since the light "flavor nonsinglet" quark does not couple to the heavy 
flavor, we see q1 ns (n, Q 2 , -P 2 )|u g ht = QL n s( n ' Q 2 >P 2 )- When all the flavor quarks are light, we usually treat the quark 
sector which consists of the flavor-singlet and nonsinglet combinations defined as follows: 

n f n f 7 

<£ = £«*> & fl ^ ( A1 ) 

i=l i=l ' 

Also we introduce the following quark-charge factors: 

1 «/ , "/ 

( e2 ) = -£ e ?> (e 4 ) = -E e '- (A2) 
nf ■* nf 

1 i=l 1 i=l 

These parton distributions, q^, qJj S and G 7 |n g ht, have been investigated in Rcfs. [12, IH, [43j]. The quark parton 
distributions qlJiight, <Z#| light and q\ ns are related to q~g and qJj S - Indeed they are expressed in terms of q$ and qJj S 
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n f =4, Q=30GeV , P =0.35GeV , A=0.2GeV, m c =1 .3GeV, MS 
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n f =4, Q 2 =30GeV 2 P 2 =0.35GeV 2 A=0.2GeV, m.=1.3GeV, MS 
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FIG. 4: Parton distributions in the photon in MS scheme for Q 2 = 30GeV 2 . The other parameters are the same as in Fig. [3] 



as follows: 



2 D 2\ 



q]i{n,Q ,P )L 



light 



I r .2 



q}(n, Q\ P 2 ) + *l f 2 \ 2 ql s (n, Q 2 , P 2 ) 



(e 4 ) - {e 2 Y 



ll ns (n,Q 2 ,P 2 



nj - 1 (e 4 ) L - (e 2 ) 



L n l 



ql s (n,Q',P'), 



(A3) 
(A4) 

(A5) 



where the NLO expressions of qg(n, Q 2 ,P 2 ) and qJj S (n,Q 2 , P 2 ) are given in Eqs. (2.33), (2.35) and (2.37) ofRef. [H]. 
The transformation rule for qj and qj^g from the MS scheme to the DIS 7 scheme are given in Eqs. (3.33) and (3.34), 
respectively, of the same reference. 
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FIG. 5: Parton distributions in the photon in DIS 7 scheme for nj = 4, Q 2 — 5GeV 2 , P 2 = 0.35GeV 2 with m c = 1.3GeV and 
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FIG. 6: Parton distributions in the photon in DIS 7 scheme for Q 2 = 30GeV 2 . The other parameters are the same as in Fig. [S] 
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